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We show that high resolution Resonant Inelastic X-ray Scattering (RIXS) provides direct, element- 
specific and momentum-resolved information on the electron-phonon (e-p) coupling strength. Our 
theoretical analysis demonstrates that the e-p coupling can be extracted from RIXS spectra by 
determining the differential phonon scattering cross section. An alternative, very direct manner 
to extract the coupling is to use the one and two-phonon loss ratio, which is governed by the e-p 
coupling strength and the core-hole life-time. This allows measurement of the e-p coupling on an 
absolute energy scale. 



PACS numbers: 78.70.Ck 

Often novel electronic properties of a material can 
be understood by systematically unravelling the interac- 
tion between its electrons and phonons. Tunable elec- 
tric transport properties in molecular crystals, for in- 
stance, are explained by the presence of a strong electron- 
phonon (e-p) coupling pQ. The dressing of electrons by 
phonons is also responsible for the colossal magnetoresis- 
tance effect in manganites [2]. More delicate is the role 
that the e-p interaction plays in high T c superconducting 
cuprates - topic of a persisting debate [3H5] . The lack of 
a technique to measure the e-p coupling strength perpet- 
uates this controversy. Here we show that high resolution 
Resonant Inelastic X-ray Scattering (RIXS) can fill this 
void as it gives direct, element-specific and momentum- 
resolved information on the coupling between electrons 
and phonons. We provide the theoretical framework re- 
quired to distill e-p interaction strengths from RIXS, par- 
ticularly in strongly correlated transition metal oxides 
such as the High T c cuprates. 

In RIXS experiments one scatters high energy, x-ray 
photons inelastically off a material [6]. The energy of 
the incident photons is chosen such that it coincides, and 
thus resonates, with an intrinsic electronic excitation of 
the material under study - one of the materials' x-ray 
absorption edges. At present the highest energy resolu- 
tions are reached at the L-edge of transition metal ox- 
ides, where an incident photon launches a 2p electron 
out of the atomic core into an empty 3d state around the 
Fermi-level. This highly unstable intermediate state de- 
cays rapidly, typically within 1-2 femtoseconds, so that 
the 2p core-state is refilled and an outgoing photon emit- 
ted. The state-of-the-art resolution is such that photon 
energy loss features on an energy scale of 25 meV can be 
distinguished at a copper or nickel L3-edge [7HS]- 

This resolution has brought phonons within the energy 
window of observation and indeed last year for the first 
time phonon loss features were resolved in RIXS (BlU0|. 



To put this achievement in perspective, one should real- 
ize that the incident photons at the Cu L-edge have an 
energy of around 930 eV, implying experiments have a 
resolving power better than 10 4 . Advanced instrumen- 
tation will drive this up further. Here we show how 
the progress in accuracy allows the extraction of a num- 
ber of characteristics of the e-p interaction directly from 
RIXS, including spatial information on the e-p coupling 
strength. No other experimental technique has access to 
such e-p characteristics, particularly in strongly corre- 
lated 3c? transition metal oxides. The advantage of RIXS 
is that x-ray photons carry an appreciable momentum. 
During the scattering process they can transfer it to a 
phonon in the solid, so that RIXS can sample the phonon 
dispersion. Present methods to measure the e-p inter- 
action do not have any access to q-dependent informa- 
tion. Electron tunneling for instance, lacks momentum- 
dependence and thus cannot probe the strength of the 
electron-phonon coupling at a given wave vector. More- 
over, it is intrinsically surface sensitive and suffers from 
the practical difficulty to make good yet partially trans- 
parent barriers [11] . Another asset of RIXS is its element- 
specificity: in copper-oxides phonons can not only be ac- 
cessed at the Cu L-edge, but also at the O if -edge - 
resolution permitting. 

Measuring phonon dispersions is an interesting new 
utilization of RIXS but as a technique it is up against 
inelastic neutron or (nonresonant) x-ray scattering. We 
make the case here, however, that what really sets RIXS 
apart is its capability to measure momentum dependent 
e-p couplings. To understand this capacity we start 
by introducing the e-p coupling Hamiltonian in generic 
form [12] 

H°~* = ^ x dt-M b U + b -^)- (!) 
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The Hamiltonian couples with a strength M q A phonons 
with momentum q and branch index A, which are created 
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by the operator 6 qA , to electrons that are described by the 
fermionic operators d^. The phonon dynamics are gov- 
erned by the phonon Hamiltonian H p = E q A u^b^b^. 

For definiteness we focus on the cuprates which have 
all their Cu 3d orbitals filled apart from a single hole 
in the x 2 -y 2 orbitals, but the analysis presented here is 
more general [T3]. In the Cu L-edge RIXS process, the 
x 2 -y 2 state is transiently occupied as an electron from 
the core is launched into it. As mentioned above, the 
filled x 2 -y 2 state is very short lived and quickly decays 
again. There is then a certain probability that after the 
decay process a phonon is left behind in the final state. 
Obviously this probability is related to the coupling of 
the x 2 -y 2 electron to this particular phonon. Note that 
the presence of the core-hole ensures that the interme- 
diate state is locally charge neutral. The phonons that 
couple to the 3d x 2 -y 2 quadrupolar charge moment ap- 
pear as loss-satellites in L-edge RIXS, due to the absence 
of a charge monopole. To understand how the electron- 
phonon coupling is precisely reflected in the RIXS inten- 
sity we evaluate the RIXS amplitude for phonon scat- 
tering Aq from the Kramers-Heisenberg equation [14] 
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where Wd t is the detuning en- 
ergy of the incident photons from the L-edge resonance, 
and E n is the energy of intermediate state n with re- 
spect to the resonance. The initial, intermediate and 
final states are \i), \n) and |/) and the momentum loss of 
the photons in the scattering process is q. The dipole op- 
erators D ~ p • A in the scattering amplitude first create 
a 3d electron out of the localized 2p core-state and then 
annihilate it. Without polarization details we can denote 
it as D = y R {e- lki " R d R p R + e ik °» tR d R p R ). An impor- 
tant quantity is the inverse core-hole life-time T, which 
determines the ultra-fast 1-2 fs time scale for the RIXS 
process. 

It is well-known that due to the presence of the in- 
termediate states, it is in general impossible to evaluate 
RIXS scattering intensities exactly, even in model sys- 
tems. One therefore often resorts to finite-size cluster 
calculations to compute RIXS spectra [13 [H] • The e-p 
scattering problem at hand, however, is an exception. It 
can be solved exactly by means of canonical transforma- 
tions. The fact that we are dealing with harmonic bosons 
allows for such a solution. We start by considering the 
limit of the coupling to an optical Einstein phonon and 
subsequently generalize to multiple dispersive phonons. 

Einstein phonon For a local, non-dispersive Ein- 
stein phonon the e-p Hamiltonian reduces to H = 
E R Mdld R (bl + b R ) + Er^rV- This Hamilto- 
nian can be diagonalized by a canonical transforma- 
tion [12] H = e s He- s 1 with 5 = ^ R S R where Sr = 
^r^r^(^r — ^r)- As t nere is a single core hole present, 



this results in H 
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where the last 



term merely shifts the energy at which the resonance 
occurs. The phonon eigenstates are all the possible 
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FIG. 1. Resonant Inelastic x-ray Scattering amplitude for 
phonon loss for Einstein phonon with energy uio for different 
values of the dimensionless electron-phonon coupling M/uiq. 
The inverse core-hole lifetime of T/too = 5 is a typical value 
at the copper L-edge. The RIXS amplitude is evaluated at 
the incident energy that corresponds to the maximum in the 
x-ray absorption signal. 



occupations {na(m)} of the phonon states with ener- 
gies E m = J2 R n R (m)ujo — M 2 /uq. The eigenstates 
|V>m.) of H are related to the eigenstates \ip m ) of H by 
\^m) — |f/> m ). Inserting this transformation in the 
scattering amplitude Aq gives 
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with 



Wdot + iT, 



and n R , 



71r, and n' R are the oc- 



cupations in the ground, intermediate, and final states, 
respectively. The transferred momentum q = ki n — k out 
springs from the dipole operators. The Franck-Condon 
overlap factors in the numerator can be evaluated ex- 
actly and are conveniently expressed in the functions 

B mn (g) = (-lrV^BElU [n [ :$^ZZi)v 

g = M 2 /luq. The RIXS amplitude to excite n' phonons 
when starting from the ground state is 
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The resulting RIXS spectra for a typical weak, interme- 
diate and strong coupling case are shown in Figs, [l] and [2] 
It is clear that for stronger e-p interactions, a larger num- 
ber of multi-phonon satellites carry appreciable weight. 
However, the first important observation is that the am- 
plitude of the zero-loss peak A^ is absolutely dominant 
(see Fig. [2]) . This is even so in the strong e-p coupling 
regime, as long as M/Y < 1. As at the Cu L3-edge 
r = 280 meV [IT], this corresponds to the physical situ- 
ation. This observation is empirically supported by the 
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(a) M/6>„=4, I7m =0.5 (b) M/w =4, r/m =5 




(c)MM=UM 




FIG. 2. Calculated RIXS amplitude for phonon loss as 
a function of loss energy lui OS s and incident energy Wj n = 
tt'dei + M 2 /ujo in the case of (a) strong coupling and very long 
core-hole lifetime, (b) strong e-p coupling, and (c) intermedi- 
ate/weak e-p coupling. 

fact that a dispersion of magnon [5] and bimagnon 
excitations have been observed in L- and K-edge RIXS. 
Such dispersion cannot be present if these magnetic ex- 
citations are always accompanied by (multiple) phonon 
excitations that carry away momentum. 

The exact amplitude for exciting a single phonon in the 
RIXS process is A« = (e~ a / 1 yfg)Y^ =0 g n {n- g) /[n\{z + 
(g — n)u)o)]. In leading order in the e-p coupling constant 
this amplitude is A^ = M/z 2 . The differential single- 
phonon RIXS scattering amplitude is therefore di- 
rectly proportional to the e-p coupling constant M, see 
Fig. [3] By increasing z and thus tuning away from the 
absorption edge the scattering amplitude is reduced. 

Dispersive phonons. In close analogy with Eq. (J2J), the 
case of multiple dispersive phonon modes can also be 
solved exactly. The scattering amplitude is 

^ = E e<qR E 

R m L 
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where n kA , n kA , and n kA are the occupation numbers of 
the modes indexed by k and A in the ground, intermedi- 
ate, and final states, respectively. The sum over m is over 
all intermediate state occupations. Together with the 



Multi-phonon RIXS ratios 
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FIG. 3. Relation of the multiple phonon resonant RIXS am- 
plitude to the electron-phonon coupling strength M/r. The 
ratio of the n + 1 phonon yl( n+1 ' to n phonon loss am- 
plitudes is shown. F is the inverse core-hole life-time. In the 
physical relevant regime these curves do not depend on the 
phonon frequency ujq. The UCL expansion (straight gray line) 
gives accurate results for M/Y < 0.2. 

Franck-Condon factors this constitutes an exact, closed 
expression for the RIXS response. 

UCL approximation Even when the e-p coupling g is 
not small, one can obtain the RIXS amplitude in ap- 
proximate form by using the fact that the time-scale of 
a typical phonon (80 meV ~ 52 fs) is much slower than 
the ultrashort RIXS time-scale (1.6 eV ~ 2.6 fs at the Cu 
K-edge). This separation of time-scales suggests that the 
scattering process contains a viable expansion parameter, 
which is small even for a fast phonon. The Ultra-short 
Core-hole Lifetime (UCL) expansion formalizes this ob- 
servation [3TJ [52] . In the present case it boils down to an 
expansion of the Kramers-Heisenberg expression in terms 
of M qA /r. With the UCL expansion we obtain a com- 
pact, approximate expression for the phonon scattering 
amplitude that is also valid at finite temperature. In the 
case of dispersive phonons, Eq. Q is expanded as 

- oo 

>W i «-X> <q " R II E «x\e- SR ^\n kX ) 

R k,A n kA =0 

oo 

x (n kA | e s ^ |n kA ) E (Ka - gkx)wk\/z) 1 (5) 

1=0 

where S RkA = 4^R^fe lk R (6 t _ kA - 6 kA ) and 5kA - 
|M kA /o; kA | 2 . As z is large compared to the phonon and 
e-p energy scales, we retain the leading terms in I and 
find the inelastic RIXS amplitude 

^hcL = ^(f\f>U + b_ cl>x \i). (6) 

to first order in M kA /z. Again, multi-phonon contri- 
butions to Aq are suppressed proportional to M kA /T. 
This result implies that momentum dependent RIXS can 
directly map out the q-dependence of the e-p coupling 
strength M qA . From this information for instance the 
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spatial range of the e-p interaction M r \ can be deter- 
mined as it is directly related to the Fourier transform of 

MqA- 

RIXS provides still another method to extract e-p 
coupling strengths, which is expected to be particu- 
larly powerful in the case of weakly dispersive optical 
phonons, for instance modes around 80 meV in the high 
T c cuprates [23j EI]- From the UCL expansion one 
finds that the ratio of the one and two-phonon loss am- 
plitude also directly reflects the e-p coupling constant: 
AW/AW = \J%Mjz. This is confirmed by the exact so- 
lution, for M/T <C 1, see Fig. [3} For strongly dispersive 
phonons this method also works if both Aq^ and A^ are 
measured throughout the Brillouin zone. 

In the analysis above we concentrated on transition 
metal L-edge RIXS, which has the advantage of a photo- 
electron launched directly into the 3d state. A certain 
3d orbital can be selected by choosing the polarization of 
incident and outgoing x-rays [23 [26J , so that e-p char- 
acteristics related this particular 3c? orbital can be mea- 
sured [5]. It is straightforward to include d-d excitations 
that couple to breathing and Jahn- Teller phonons [57], 
which is of great importance in the study of Jahn- Teller 
polarons. Also at the O K-edge phonons can be probed 
but as this edge is at lower energy, the photons have less 
momentum and a smaller part of the Brillouin zone can 
be probed, which also holds for Cu M-edges. As hard 
x-ray transition metal K-edges do not suffer this disad- 
vantage they do provide in principle a viable method to 
measure extensively momentum dependent phonon prop- 
erties and e-p interactions [THl EB] • 

The framework presented here to extract the e-p cou- 
pling interaction bestows on RIXS the unique poten- 
tial to provide direct, element-specific and momentum- 
resolved information on the interaction between electrons 
and phonons on an absolute scale. In weakly correlated 
electron systems these properties can be computed with 
modern ab initio electronic structure methods, for in- 
stance in the newly discovered iron pnictide superconduc- 
tors [22]; and our framework to distill them from RIXS 
allows a direct comparison. In strongly correlated mate- 
rials, particularly the high T c cuprates, these assets make 
high resolution RIXS a unique tool to unravel the inter- 
action between its electrons and phonons. 
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